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Phonon lifetime calculations from first principles usually rely on time consuming molecular dy- 
namics calculations, or density functional perturbation theory (DFPT) where the zero temperature 
crystal structure is assumed to be dynamically stable. Here is presented a new and effective method 
for calculating phonon lifetimes from first principles. This method is not limited to crystallographic 
phases stable at K, and provides a scheme more effective than most corresponding molecular dy- 
namics calculations. The method is based on the recently developed self consistent ab initio lattice 
dynamical method (SCAILD) and is here tested by calculating the bcc phase phonon lifetimes of 
Li, Na, Ti and Zr, as representative examples. 



Calculations of phonon lifetimes are of great interest 
to the materials science community, partly in that it al- 
lows for the prediction of thermal transport properties, 
which in themselves are important from an engineering 
perspective, but also, since these properties are closely 
interrelated to perhaps the even more intriguing thermo- 
electric effect pp. 

Historically phonon lifetime calculations from first 
principles have been performed for quite some time, for 
instance see the work of Katsnelson et al. 2\ or De- 
bernardi et al. [3J. More recently Koker [3] successfully 
calculated the thermal transport coefficients for MgO 
within the context of ab initio molecular dynamics (ab- 
MD) [5], and Bonini et al. estimated the optical phonon 
line-widths in graphene using density functional pertur- 
bation theory (DFPT) [HHE]- 

Up to present date ab-MD and DFPT have been the 
two main methods of choice used for calculating phonon 
lifetimes from first principles. However, the use of ab-MD 
suffers from the long simulation times required to sam- 
ple the correlation functions, whereas DFPT is limited 
to crystal structures stable at zero temperature, in that 
it only calculates third order anharmonic contributions. 
Thus, there is still a need for more effective methods, 
preferably based upon schemes that also allow lifetime 
calculations in crystallographic phases stabilized only at 
finite temperatures, such as the bcc phase of Ti, Zr and 
Hf and the cubic phases of numerous shape memory al- 
loys. 

In this paper a novel and effective method for calculat- 
ing phonon-lifetimes, which do not suffer from the lim- 
itations of ab-MD and DFPT, will be presented. The 
method is based on the recently developed self con- 
sistent ab initio lattice dynamical scheme (SCALD) [9], 
which has previously successfully been used in predict- 
ing phonon frequency renormalization by strong anhar- 
monicity [THl E] • The newly developed method will here 
be used to calculate phonon lifetimes in the bcc phase of 
Li, Na, Ti and Zr, in order to illustrate its effectiveness. 

The SCAILD method is based on the calculation of 
Hellman-Feynman forces of atoms in a supercell. The 
method can be viewed as an extension of the frozen 



phonon method |12j . in which all phonons with wave 
vectors q commensurate with the supercell are excited 
together in the same cell by displacing atoms situated at 
the undistorted positions R + b^, according to R + b^ — > 
R + bo- + Ur ct , where the displacements are given by 

u ^ = -/^E^<s eiq(R+bCT) - C 1 ) 

Here R represent the N Bravais lattice sites of the su- 
percell, b CT the position of atom a relative to the lattice 
site, e qs are the phonon eigenvectors corresponding to the 
phonon mode, s, and the mode amplitude is calcu- 
lated from the different phonon frequencies w qs through 

where n qs = n(^- 3 ^), with n(x) = l/(e x — 1), are the 
phonon occupational numbers, M a the atomic masses 
and T is the temperature of the system. The phonon 
frequencies, u; qs , are obtained through the projections 
of the Fourier transformed atomic forces, F q , onto the 
eigenvectors of the corresponding mode 

Due to the simultaneous presence of all the commen- 
surate phonons in the same force calculation, the inter- 
action between different lattice vibrations are taken into 
account and the phonon frequencies given by Eq. ^ are 
thus renormalized by the very same interaction. 

By alternating between calculating the forces on the 
displaced atoms and calculating new phonon frequencies 
and new displacements through Eqs. ([l])-((3| the phonon 
frequencies are calculated in a self consistent manner. 

During the course of a SCAILD calculation frequency 
distributions, 2?(o; qs ), are generated for each of the 
modes having q-vectors commensurate with the super- 
cell [Ml [14]. By calculating the first moment, f2 qs , of 
these distributions the renormalized phonon frequencies 
are obtained [3]. 
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The novel implementation of the SCAILD scheme uti- 
lizes the possibility of extracting the phonon lifetimes 
from the frequency distributions, X>(u; qiS ). To estab- 
lish the connection between the phonon lifetime, r qs , 
of a mode and its corresponding distribution T>(cj^ s ), 
we start from the newtonian equations of motion of a 
stochastically damped phonon mode 



dt 2 ~ qs dt 



- O A" 



a(t), 



(4) 



where r qs = l/r qs is the phonon linewidth (damping), 
and a(t) is the stochastic acceleration arising from the 
interaction between the mode and the other phonons of 
the system. The replacement of interaction terms by 
stochastic variables, is in path- integral approaches as- 
sociated with the Stratonovich-Hubbard transformation 
|15j . and in more classical contexts, associated with dis- 
cussions concerning the fluctuation dissipation theorem 

Following a path similar to what is generally used in 
proving the fluctuation dissipation theorem [161 117] , the 
damping r qs is formally connected to the acceleration 
correlation function, C(t) = (a(to)a(to +t)), through the 
integral equation 
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c*(r q 



(5) 



where the brackets (, ) denote an ensemble average, Uq S 
O qs - (r qs /2) 2 and C* is given by 



dt 



C*(r qs ) = 

cos(i/ qa t) - a qs cos(i/ qs |i| + S v ) C(t)e~ r « M .(6) 



Here tan(5q S ) = 2z/ qs /r qs and a qs = (r qs /2)fi qs 1 . 

Even though Eq. (5][6) in principle provides a connec- 
tion between the stochastic friction forces and the phonon 
lifetimes, we will here instead follow an alternative ap- 
proach in order to extract the phonon lifetimes, which 
is more convenient when employed in the context of the 
SCAILD scheme. 

The phonon frequencies obtained through the SCAILD 
scheme, Eq. ([3]), can be related to the homogeneous 
(transient) solutions to Eq. Q, i.e solutions achieved 
with a(t) — 0, through 



1 d 2 A° 



»/4 rt „ dt «/4 „ „ dt 



r d A a 
i qs "-^i qs 



fit 



(7) 



qs 



By assuming that the temperature is high enough for a 
classical description to be valid, we can use the transient 
solutions of Eq. Q together with Eq. (JtJ) to obtain the 
following relation 



dt 



(r. 



) 2 + « s 



(8) 
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FIG. 1: (Color online) Calculated phonon lifetimes (full black 
circles) for (a) bcc-Li and (b) bcc-Na, displayed together 
with experimental data (red squares) |21H23|. The full lines 
connect lifetimes corresponding to longitudinal modes, the 
dashed lines connect lifetimes corresponding to the transverse 
Ti modes and the dotted line connects lifetimes corresponding 
to the transverse T2 modes. The calculations were performed 
at a temperature of 293 K , whereas the experimental data 
for bcc-Li and bcc-Na were obtained at 293 K and 296 K, 
respectively. 



between the frequency distributions and the phonon life- 
times. Thus it becomes evident that the distribution 
^ calculated within the SCAILD scheme corresponds, 
within a multiplicative constant, to the dynamical struc- 
ture factor of a stochastically damped harmonic oscillator 
in the high temperature classical limit [HI [19] . 

Another important observation to make at this point, 
is that the frequency distributions can be calculated de- 
spite the absence of any explicit time dependence in the 
SCAILD scheme, by virtue of the stochastic damping and 
the ergodic principle. 

Regarding the other computational details of the force 
calculation the VASP package [20] was used, within the 
generalized gradient approximation (GGA). The PAW 
potentials required energy cutoffs of 210 eV, 125 eV, 200 
eV, 170 eV for Li, Na, Ti and Zr, respectively. The k- 
point mesh used was 6x6x6 Monkhorst-Pack, together 
with a Methfessel Paxton smearing of 0.2 eV. The super- 
cclls used were obtained by increasing the bec primitive 
cells 8 times along the primitive lattice vectors, resulting 
in 512 atom cells. 

The phonon frequency distributions were obtained 
through 200 SCAILD iteration together with a gaussian 
smearing of 0.05 THz 2 applied to each of the sampled 
squared frequencies. 
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FIG. 3: (Color online) Calculated phonon lifetimes (full black 
circles) for (a) bcc-Ti and (b) bcc-Zr, displayed together with 
experimental data (red squares, blue diamonds and green 
triangles) 24, 25 . The full lines connect lifetimes correspond- 
ing to longitudinal modes, the dashed lines connect lifetimes 
corresponding to the transverse Ti modes and the dotted line 
connects lifetimes corresponding to the transverse T2 modes. 
The calculations were performed at a temperature of 1300 
K, whereas the temperature for which the experimental data 
were obtained is given in the figure. All the experimental 
lifetimes correspond to transverse Ti modes. 



FIG. 2: (Color online) The calculated Ti mode frequency 
distribution for bcc-Li and bcc-Zr along the [££0] direction, 
for £ = 0.125,0.250,0.375,0.500 (dashed black lines). The 
full red curves were obtained through the fitting of the data 
with Eq. Q. The calculations were performed at the finite 
temperature of 293 K. 



In Fig. [2] and Fig. [3] the calculated frequency distri- 
butions of bcc-Li and bcc-Ti are presented together with 
their respective fits to Eq. ffify. Here the distributions 
obtained from the SCAILD iterations are found to be in 
excellent agreement with the fit to Eq. pL Furthermore, 
in Fig. [I] and Fig.[3]the phonon lifetimes extracted from 
the fits to Eq. (|8| are displayed together with experi- 
mental data. Here good agreement between theory and 
experiment can be found in the [110] direction for bcc- 
Na, bcc-Ti and bcc-Zr, whereas a potentially quite huge 
discrepancy can be found for bcc-Zr in the [112] direction 
at £ = 0.35. This discrepancy is referred to as potential 
due to the fact that the experimental q-point correspond- 
ing to £ = 0.35 does not belong to the q-point set used in 
the calculation, and because linear interpolation between 



lifetimes corresponding to different q-points is not always 
a good approximation to employ. Finally, it can also be 
observed that the theoretical phonon lifetimes in the case 
of bcc-Li are well within the experimental error bars. 

In summary, a novel, effective and simple method 
for calculating phonon lifetimes have been presented 
together with some test calculations on the bcc phase 
occurring at room temperature as well as for tempera- 
tures well above 1000 K, for Li , Na, Ti and Zr. The 
method has provided good estimates to most of the 
phonon lifetimes, and show promise as an effective 
tool in calculating phonon line-widths, especially for 
crystal structures that can only be stabilized at ele- 
vated temperatures. Thus opening the door to phonon 
lifetime calculations from first principles in cases when 
anharmonic interactions of order > 3 are important for 
achieving a correct description of the system under study. 

I would like to thank the Swedish National Infrastruc- 
ture for Computing (SNIC) for making this work possible 
through the allocation of computational time to account 
SNIC 006/10-3 at the PDC Lindgren cluster. 
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FIG. 4: (Color online) The calculated Ti mode frequency 
distribution for bcc-Ti and bcc-Zr along the [££0] direction, 
for £ = 0.125,0.250,0.375,0.500 (dashed black lines). The 
full red curves were obtained through the fitting of the data 
with Eq. Q. The calculations were performed at the finite 
temperature of 1300 K. 
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